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The ambiguity of the structure of !?Be especially in the configuration of !?Be ground state has attracted a lot 
of attention recently. We notice that the nuclear reaction cross section oR at low energy region is sensitive to the 
surface structure of !?Be, which is greatly impacted by the ground state configuration of !?Be especially by the 
occupancy probability of the s orbital component. By using existed interaction cross section data of !?Be on C 
at 790 MeV/nucleon and Glauber model, the upper limit of the s orbital occupation probability of '?Be ground 
state is roughly determined to be about 5696 with Single Particle Model calculations. This demonstrates that 
the method is very promising to determine the s orbital component of !?Be with proper nuclear-matter density 
distribution calculations for different orbitals of !?Be ground state. Hence we bring forward to determine the 
s orbital component of !?Be by measuring the og of !?Be on C and Al at several tens of MeV/nucleon. In 
this paper, the feasibility and detailed experimental scheme of the or measurement are carefully studied. The 
precision of the s orbital occupation probability of ‘Be ground state is expected to achieve 9% by using the 


proposed 2% op data. 
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I. INTRODUCTION 


During the past decades, studies on exotic nuclei have been 
stimulated considerably owing to the enormous development 
of radioactive ion beam (RIB) technique. Many peculiarities 
of exotic nuclei have been revealed, such as the halo/skin-like 
structure [1], the cluster structure [2] and the breakdown of 
shell closure [3, 4]. Studies on these peculiarities greatly im- 
proved understanding of the exotic nuclei, e.g. the beryllium 
isotopes. Theoretical researches on properties of Be isotopes 
are based on the three-body model [5, 6] and the density- 
dependent relativistic mean-field model [7, 8]. The neutron 
halo of !^Be has been well explained using these models [5— 
8]. Located between the halo nuclei !! Be and ‘Be, !?Be is 
an interesting combination of the peculiarities and plays a key 
role in the beryllium chain. In shell model, !?Be is supposed 
to be a “magic nucleus" with simple structure. But recent ex- 
periments [4, 9, 10] provided direct evidence for breakdown 
of the N = 8 shell closure in !?Be and the (s, d) intruder 
states. In principle this intruder-state configuration can cause 
halo-like structure, yet neither the wide momentum distribu- 
tion [11] nor the relative large two-neutron separation energy 
(Son = 3.67 MeV [12]) indicates the signature for a halo. 
There exists controversy between the recent theoretical and 
experimental results on halo-like structure of !?Be. The gi- 
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ant deformation in the ground state of !?Be [13] is also pre- 
dicted by Antisymmetrized Molecular Dynamics (AMD), the 
ground state structure of !?Be is still ambiguous. Also, ex- 
plorations in the level structure, parity, spin and deformation 
of the excited states of !?Be with molecular cluster model 
draws quite a lot of attentions [14-18]. The ground state 
structure properties of !?Be are indispensable in the studies 
because it provides fundamental information for the studies 
of the excited states. Therefore, it is of significance to study 
the ground state structure of !?Be. 

As an essential property of the ground state structure, the 
nuclear-matter density distribution of !?Be not only provides 
basic structure information such as nuclear-matter radius, but 
also helps to determine whether the !?Be ground state has a 
halo structure. Recently, particle-particle random-phase ap- 
proximation (pp-RPA) [19] and microscopic no-core shell- 
model (NCSM) calculations [18] have shown that the ground 
state wave functions of !?Be are dominated by the p shell 
configuration, which is in conflict with the previous calcu- 
lations [20—24] and the knockout measurements [4, 9]. The 
nuclear-matter density distribution may be used to determine 
the configuration mixing between the 1p; 5, 1d5/2 and 25; /2 
orbitals, thereby help to resolve the inconsistencies between 
current data and various theoretical models. Therefore, we 
were motivated to determine the nuclear-matter density dis- 
tribution of !?Be ground state. 

In the following text, the experimental method of the opr 
measurement is introduced in Sec.II, the feasibility analysis 
and optimization of the reaction cross section (o R) measure- 
ment are elaborated in Sec.IIL, the way to determine the s 
orbital component of !?Be ground state through op is illus- 
trated in Sec.IV, and a summary is given in Sec.V. 
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Il. GENERAL SCHEME 


Typically the nuclear-matter density distribution is investi- 
gated by measuring the øp (or the total interaction cross sec- 
tions a). This involves theoretical models, and several meth- 
ods were developed to study the total reaction cross section, 
such as the multi-step scattering theory of Glauber [25], the 
transport model method of Ma et al. [26, 27] and the semi- 
empirical formula of Kox et al. [28] and Shen et al. [29] 
A series of investigations to !?Be were carried out in the 
past decades. By measuring the o; on Be, C and Al at 
790 MeV/nucleon, Tanihata et al. successfully determined 
the effective root-mean-square (RMS) radius of !?Be through 
a Glauber model in 1988 [30]. Liatard et al. measured 
the og of radioactive Be isotopes on Cu at around 25 to 
65 MeV/nucleon and deduced the radii by using a simple mi- 
croscopic model [31]. Later, Warner et al. measured the og 
of !?Be on Pb and Si at about 30 to 60 MeV/nucleon and ob- 
tained the radius of !?Be [32]. However, the studies have 
several drawbacks. First, their energy region did not cover 
both low and high energy regions. Data at low energy region 
give constraint to the outer structure, and data at high energy 
region provide more information on the core part, hence low 
and high energy region data are both needed to obtain accu- 
rate nuclear-matter density distribution. Second, generally, 
their targets were too heavy. This makes øp insensitive to the 
surface structure of !?Be. Third, their models were too simple 
to interpret the data especially at low energies. The results did 
not give detailed nuclear-matter density distribution of !?Be 
especially for the surface area. 

In obtaining detailed outer structure of the nuclei, two 
improvements were made recently for extracting accurate 
nuclear-matter density distribution of nuclei. First, the 
method of proton elastic scattering at intermediate energies 
was developed. Ilieva et al. applied proton-scattering method 
in 2012 to !?Be [33] and determined the nuclear-matter den- 
sity distribution of !?Be ground state. But the result has con- 
siderable uncertainty because of various parametrizations. So 
the method is not well established yet for studing the surface 
structure of especially unstable nuclei. Second, the applica- 
bility of Glauber model in the whole energy region was stud- 
ied and a Modified Optical Limit Glauber model (MOL[FM ]) 
was developed by incorporating the Fermi motion of nucle- 
ons in the finite-range MOL [34]. The MOL[FM] has re- 
duced the discrepancy between the data of calculation and 
measurement to just 176—246 for the whole energy region. By 
measuring the og on Be, C and Al targets at intermediate 
energies, the nuclear-matter density distributions of ??C [35] 
and !"Ne [36] were determined accurately with MOL[FM]. 
Through this method we obtained the nuclear-matter density 
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distribution of *Li [37] precisely and the relevant paper about 
our further study is in preparation. Therefore, MOL[FM] 
provides a powerful tool for interpreting the o; or op data, 
and measuring the c; or og determining the nuclear-matter 
density distribution with MOL[FM] is still a good method at 
present. 

Thus, we were motivated to precisely measure op of !?Be 
at low energy region and extract the nuclear-matter density 
distribution of !?Be ground state using MOL[FM]. 

Refer to the measurements of cg, transmission method is 
usually used. Typical experimental procedures are given in 
Ref. [38]. The op is obtained by Eq. (1), 


on(B) - -, nR, () 


where F is the energy point; t is the target thickness ex- 
pressed by the target particles numbers per unit area; and 
R = Nou/Nin is the ratio of outgoing projectile particles 
number to incident projectile particles. Since energy of the 
projectile particles decreases as passing through the target, 


we determine the energy point of op by mean energy Emean, 
which is given by 


t 
E 
mean — | Un dz, 
0 t 


(2) 
where f is the target thickness, and E(x) is the residual en- 
ergy of incident beam travelling along the path by distance 
x. E(x) can be calculated by the improved Bethe-Bloch for- 
mula [39]. 

The main relative error of ap can be written as 


Aog M, | AtY | pt ARY | AR\* 
OR Xt In R R in FE 
(3) 

where At/t stands for uncertainty of the target thickness, the 


subscripts “sys” and "stat" denote the systematic and statis- 
tical error of R, respectively, and the statistical error of R is 


given by 
( ) tat 
R S 


because it follows the binomial distribution. In practice, 
Rin/ Rout is taken as R in order to remove the events inter- 
act outside the target, Ri, and Rou are the ratios of Nou/ Nin 
corresponding to the target-in and target-out measurements, 
respectively. Accordingly, the relative error of ap becomes 


1-R 
Nin’ 


(4) 


r= Rin 1— Rout 
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Besides the contribution of statistical error, the systemic 
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uncertainty of or is mainly from the correction of the num- 
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ber of inelastic-scattering events merged into the non-reaction 
events. By using the method in Ref. [40] and Monte Carlo 
simulation, the systematic error of c a can be limited within 
196—296 [34]. So 196—296 total uncertainty can be achieved 
1f sufficient reaction events are recorded. In fact, this was 
achieved recently in most experiments of the kind. 
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Fig. 1. (Color online) Current results of the nuclear-matter den- 
sity distributions of !?Be (XB97: [41], SI12: [33], GG05: [42], 
IT88: [30]). The solid line is given by the sum of the nuclear-matter 
density distribution of 10Be core in ref. [33] and that of two valence 
neutrons calculated in section 4 with the configuration in ref. [4]. 


Based on available nuclear-matter density distributions of 
12Be from experiments and theories (Fig. 1), we calcu- 
lated the or (Fig. 2) at different energies by MOL[FM]. 
The momentum width and the finite-range parameter 3 were 
taken from Ref. [34]. From Fig. 2 one sees that below 
~50 MeV/nucleon, cg data with uncertainty around 2% is 
sufficient to distinguish several previous results and deter- 
mine whether !?Be has a clear halo-like structure. 

However, the specific requirement on experimental condi- 
tions for a 2% precision of og is still unknown, and investi- 
gation is needed for further discussion of the method's feasi- 
bility. 


IH. FEASIBILITY STUDY AND OPTIMIZATION OF oz 
MEASUREMENT 


To make sure that the experiment is practical, we have stud- 
ied the feasibility and optimized the experimental method. 
The experiment feasibility mainly relies on detecting system, 
reaction targets and !?Be beam. In the energy region below 
50 MeV/nucleon, AE — E method is often used for particle 
identification. The required energy deposition of !?Be in the 
AE detector is about 10 MeV/nucleon, while the required en- 
ergy in the E detector shall be above 10 MeV/nucleon, just for 
ensuring validity of the particle identification. So the energy 
of !?Be right before the AF detector is supposed to be over 
20 MeV/nucleon. Usually Si detector and scintillator detector 
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Fig. 2. (Color online) The calculations of o g() corresponds to 
current nuclear-matter density distributions of 12Be (XB97: [41], 
SI12: [33], GGOS: [42], IT88: [30]). The upper group of lines corre- 
sponds to the or of !?Be + Al, and the lower group of lines corre- 
sponds to the op of '?Be + C. 


are used as the AF and E detector, respectively. Their thick- 
nesses depend on the specific energy of !?Be right before the 
AE detector. As for the particle identification before the re- 
action target, Bp-A.E-TOF technique is often used. Typical 
layout of the detecting system is shown in Fig. 3. 


<— TOF—> AE-E 


P Reaction 


FIPL F2PL target 


Fig. 3. Schematic diagram of typical experimental setup. 


Besides the detection system, the target material and thick- 
ness t and incident energy of the !?Be beam E, shall be deter- 
mined before the experiment. Regarding the target material, 
12C is a good candidate, because its nuclear-matter density 
distribution is well determined and its mass number is com- 
parable with !?Be, so its cg is more sensitive to the !?Be 
surface structure. Another target, 27 A], is also needed for re- 
ducing the target-dependence of the result. 

For the t and Ej, there are direct impact factors, such as 
the transmission rate R, !?Be outgoing energy E, after the 
reaction target, and energy point Emean with certain restriction 
for each of these factors, hence it is difficult to consider t 
and Ej, separately. We studied the relations between these 
parameters through calculations, which includes: 


(1) Fou for certain Ej and t, E, was calculated by 
LISE++ [43]; 


(2) Emean could be determined through (Fin + Eout)/2 ap- 
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proximately; 


(3) or(Emean) was calculated using MOL[FM] based on 
the nuclear-matter density distribution of the expecta- 
tion in Fig. 1; 


(4) R is calculated by Eq. (1) with certain ¢ and corre- 
sponding o g( Emean). 


Then, the trends of R varying with Ej, under different t and 
out were obtained. Fig. 4 shows the calculation results of 
1?Be4-C. Each dash line corresponds to the same t, each dash 
dot line corresponds to the same Four. Subsequently, we took 
into account the restrictions of R, Es and Emean to give 
proper t and Ein. As mentioned before, Emean should be less 
than 50 MeV/n to ensure the physical goal, and E, (viz. the 
energy of !?Be right before the A E detector) should be above 
20 MeV/n to ensure the availability of AE — E method. 

R is restricted by systematic error if the statistical error 
is small enough. It is of certain difficulty to achieve 0.1% 
systematic error for a directly measured quantity, like R. To 
ensure a < 2.5% precision of c p, the factor 1/In R before the 
systematic error of R in Eq. (5) shall be less than 25. Accord- 
ingly R should be less than 0.96. Therefore, proper ranges of 
t and E, are indicated as the hatched area in Fig. 4. 
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Fig. 4. The change trends of R with Fin under different t and Fou. 


In the energy region of 30-50 MeV/nucleon we intend to 
obtain three data points by using subtraction method. As 
shown by Conditions (1) and (2) in Fig. 5, at incident beam 
energy of E, and Eo, with target thickness of tı and t2, reac- 
tion cross sections cı and c» are measured by the transmis- 
sion method, respectively, and by adjusting E, and E5, both 
the outgoing beam energies Fouw can be of the same energy. 
Then, the reaction rate in the target in Condition (2) shall be 
equal to that of the corresponding thickness of f£» in the tar- 
get in Condition (1). By subtracting the two data, one obtains 
another c g as, 


sub in—1 
R (=) ty A = a 2), (6) 


Beam 


(2) 


Beam 


Fig. 5. (Color online) Schematic diagram of the subtraction method. 


where Rin—1 and Rin—2 denote the transmission rates of the 
target-in experiment in Conditions (1) and (2), respectively. 
The transmission rates of the target-out experiments do not 
appear in Eq. (6) because the reaction rate outside the tar- 
get is canceled between the two measurements. The error of 
this deduced op was determined in Ref. [40]. Through this 
method, three data points can be obtained with only two mea- 
surements. This greatly improves the efficiency of the exper- 
iment. The key point is that the outgoing energies Fout in the 
two measurements should be the same. And the interval be- 
tween the two adjacent incident beam energies Ej, should be 
over 15 MeV/nucleon in order to make sure the energy points 
are evenly distributed between 30—50 MeV/nucleon. 

Thus we have decided the range at Fin, and C and AI tar- 
get thicknesses, and the detecting system. Next, we are to 
determine the number of events we need according to the re- 
quirement of the statistical error, and estimate the requisite 
beam intensity. 

According to Eq. (3) the statistical error should be less than 
0.5%, so that it will not be a main error source, which gives 


1 1—-R 1 
|In R] R JN 
Since R is at most 0.96, Nin of 10° is already sufficient. So, 


it will take only 10? seconds of beam time under typical con- 
dition of 10? s^! beam intensity, which is very practical. 


« 0.596. (7) 


TABLE 1. Experimental scheme 


Beam Ei Eon Target Thickness R Nin 
(MeV/n) (MeV/n) (g/cm?) 
= 
62.7 24 C 2.33 0.86 2.9x10 
12Be Al 2.70 0.90 4.0x 10° 
42.4 24 C 0.90 0.94 6.7x 1o. 
Al 1.05 0.96 1.0x10 


To sum up, for the experiment a !?Be beam of about 
10? s~! intensity at 20-70 MeV/nucleon is needed to provide 
the projectile particles, and the outgoing energy of !?Be af- 
ter the reaction targets in the two measurements should be 
the same. In addition, a Si and a scintillator detector are 
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Fig. 6. Schematic diagram of the layout of RIBLL. 


needed for particle identification. Finally, we select two en- 
ergy points, find the corresponding target thickness, evaluate 
the R and decide corresponding requirement of Nin. The de- 
tailed experimental scheme is given in Table 1. 

We find that the Radioactive Ion Beam Line in 
Lanzhou (RIBLL) [44] is a suitable candidate for providing 
1?Be beam. The schematic diagram of the layout of RIBLL 
is shown in Fig. 6. 


The primary beam !9O** is accelerated by the Heavy Ion 
Research Facility of Lanzhou (HIRFL) and introduced to RI- 
BLL. It bombards the production target of Be at TO and gen- 
erates the secondary beam of !?Be. An AI degrader at C1 is 
used for an energy-loss analysis of the secondary beam sep- 
aration. Another Al degrader at T1 is used for energy degra- 
dation. Two slits at C1 and C2 is used for momentum accep- 
tance controlling. Two plastic scintillation counters at focal 
points T1 and T2 can provide the TOF information. Si solid 
state detector (SSD) at T2 can be used to provide the A E sig- 
nals for particle identification before the reaction target. En- 
ergy Ej of the incident beam right before the reaction target 
is determined by the Bp value of the fourth dipole magnet D4. 
Based on this configuration, the !?Be beam condition is sim- 
ulated by LISE++ using a typical 100enA primary beam of 
1808+ at 80 MeV/nucleon. Major parameters of the simula- 
tion result given in Table 2. It's worth mentioning that similar 
experiments have been performed on RIBLL since 2000 [45— 
49]. From their results we infer that typically the error of o 
is up to about 596 including 396—496 systematic error. So we 
need further consideration based on the real performance of 
RIBLL to reduce the error especially the systematic error. 

For the detecting system, typically a Si detectors of 
1500 um thickness is competent for the AE measurement, 
and a CsI(TI) scintillator detector of 30 mm thickness is ade- 
quate for the E measurement. It should be noted that the Si 
detector is made of single crystal, hence the concern of chan- 
neling effect to the particles being detected. By tilting the Si 
detector against the beam axis at a certain angle, the fraction 
of channeling events can be reduced to « 196 [40]. The tilting 
angle can be determined by studying the angle dependence of 
the channeling events in advance. 

Based on the above feasibility study, the oR measurement 
is reasonable and valid, and the requisite conditions can be 
satisfied. Through the measurement the nuclear-matter den- 


TABLE 2. Major parameters of !?Be beam condition at RIBLL 


Items* Value 1 Value 2 Value 3 
Production target thickness (mg/cm?) — 370 740 740 


D1 (Tm) 3.7594 3.5621 3.4911 
C1 degrader” thickness (mg/cm?) 810 1350 1350 
D2 (Tm) 3.5355 3.1014 3.0010 
T1 degrader^ thickness (mg/cm?) 0 540 1350 
D3 (Tm) 3.5335 2.8575 2.1478 
D4 (Tm) 3.5335 2.8575 2.1478 
Beam energy (MeV/nucleon) 62.73 42.42 24.06 
Beam rate (s71) 3940 3140 789 
Beam purity 99.1% 99.3% 98.9% 


a The CI slits are of 40 mm (H), and the C2 slits are of 40 mm x 140 mm 
(Hx V). 
> The degrader is plane-shaped. 


sity distribution of !?Be ground state can be determined and 
then used to extract the component of !?Be ground state. In 
Sec.IV, we elaborate that how to determine the s orbital com- 
ponent through op. 


IV. EXTRACTION OF THE s ORBITAL COMPONENT OF 
!2BE GROUND STATE 


As shown by the fermionic molecular dynamics (FMD) 
calculations (see the inset (b) in Fig. 10 of Ref. [33]), dif- 
ferent configuration mixings of the two valence neutrons in 
12Be lead to different nuclear-matter density distributions of 
12Be, and the difference is obviously indicated in the surface 
structure of !?Be. Inspired by this result, we bring forward a 
method to extract the ground state component of !?Be. 

In the method, we treat !?Be as a system of '°Be core 
plus two valence neutrons as usual. The !°Be determines the 
nuclear-matter density distribution of the core part of !?Be, 
while the two valence neutrons determine the outer part den- 
sity distribution. According to the intruder configuration of 
1?Be [4], the two valence neutrons are populated in 1p, /2> 
1d5,2 and 25, 2 states at certain occupation probabilities. So 
we can construct the outer structure of !?Be according to a 
certain configuration as long as a model can give preferable 
density distribution of the two valence neutrons correspond- 
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Fig. 7. Density distributions of the two valence neutrons correspond- 
ing to 1p1/2, 1d5/5 and 25, /2 states. 


ing to different states. By using an appropriate function to de- 
scribe the core structure of !?Be, we can construct a nuclear- 
matter density distribution of the ground state of !?Be. By 
adjusting the proportion of the components, we can find a 
nuclear-matter density distribution which is consistent with 
the experimental result. Then the corresponding proportion 
can provide configuration information of !?Be ground state. 

We have tried to extract the ground state component of 
12Be in order to verify the feasibility of the method. As men- 
tioned above, a precondition of the method is to obtain the 
valence neutron density distributions corresponding to 1p; 5, 
1d5,5 and 25, /2 states. Usually single particle model (SPM), 
three-body-model, cluster model, shell model, etc. are used 
to calculate the density distribution of the valence neutron. In 
this article, the SPM [35] is used. In this SPM, the Woods- 
Saxon potential, the Coulomb barrier and the centrifugal bar- 
rier are taken into account. The nuclear part of the potential 
assumed is written as 


2 u =i 
y- (-w evi gh) [1+ exp (z 23] 
(8) 


where Vo is the depth of the Woods-Saxon potential, Vj = 
17 MeV is the l- s strength taken from Ref. [50], 7;., = 1.1 fm 
is the radius of spin-orbit potential, Re = roA! (ro = 
1.2m) is the radius of the Woods-Saxon potential, and a = 
0.6 fm is the diffuseness parameter. Vo is adjusted to re- 
produce the separation energy of the valence neutron. Here 
we treat the two valence neutrons as equal and set the sep- 
aration energy of single valence neutron to be a half of the 
two-neutron separation energy of !?Be. The corresponding 
nuclear-matter density distributions of the two valence neu- 
trons are shown in Fig. 7. We can see that the 25; state 
has larger density in the surface. Although different mod- 
els will not give exactly the same density distributions of the 
two valence neutrons, we infer that our result is less model- 
dependant based on the fact that s intruder valence neutron 
configuration is the chief cause of the halo-like structure in 
light nuclei just as the cases of !! Li and !! Be. 
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Then we consider a configuration mixing of 1p /2, 1d5/2 
and 2.51 /2 for the two valence neutrons as follows: 


Pon = OP[ap, jy] + PP aas] (19 B)P[(;, 2]; (9) 


where o and £ denote the occupation probability of (1p; /2)° 
and (1d; py configuration, respectively. The integrals of the 
single configuration density and mixed density are all nor- 
malized to be two nucleons. Combining the !?Be core dis- 
tribution of Gaussian-Gaussian (GG) parametrzation given in 
Ref. [33], we calculated the corresponding op and obtained 
a range of cg accordingly, as indicated in Fig. 8. We can 
see clearly that the øp is sensitive to the s orbital occupancy, 
and it is difficult to derive the p and d orbital occupancies 
from the oR data. This is because the 25; /2 component con- 
tributes much more to the surface structure of !?Be than the 
other components. By using the c; of !?Be on C in Ref. [30], 
we extract the upper limit of the s orbital occupation proba- 
bility to be about 56%. The upper limits are indicated by the 
vertical line in Fig. 8(c). 


The result indicates that the s orbital is not a dominant 
component in !?Be ground state. It is supposed to be the 
reason why the halo-like structure in 1?Be is not evident, be- 
cause the s orbital component is the major contribution to the 
halo-like structure. Compared with the calculations given by 
Ref. [4, 20—22, 24], our result is relatively small (Table 3). Al- 
though our result is consistent with the 6 decay result given 
by Ref. [42], large uncertainty exists in both results. There- 
fore, higher precision extraction is needed. 


By comparing the increasing trend of øp in insets (c) and 
(f) of Fig. 8, one can sees that at low energy region, og is 
even more sensitive to the s orbital occupancy. This implies 
that the extraction uncertainty can be reduced if o p at low en- 
ergy region is used. So, we calculate a based on an s orbital 
occupation of 2846 (half the value of the upper limit) and a 
p orbital occupation probability of 25%, as the percentage of 
1p,/2 component was determined as 25 + 5% [10]. The re- 
sult is shown in Fig. 8(f), in which the hatched area indicates 
2% error band. Thereby we can extracted that uncertainty of 
the s orbital occupation probability is about 2396. And take 
into account that six data points will be available from the ex- 
periment introduced in previous section, the precision of the 
extracted s orbital occupation probability is expected to be 
around 9%. This makes the experimental measurement of o. 
at low energy region more desirable. 


Since the percentage of the 1p; /2 component has been de- 
termined by Gamow-Teller transition strengths [10], if either 
1ds,2 or 28/2 component is extracted precisely, the ground 
state configuration of !?Be shall be reliably determined. Our 
method provides a promising approach to extract the 25, /; 
component of !?Be ground state and to determine the ground 
state configuration of !?Be. Through the proposed c; mea- 
surements, we expect to extract the s orbital occupation prob- 
ability with 9% uncertainty. 
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Fig. 8. The reaction cross section (cg) of !? Be calculated at 33.2 MeV/nucleon and 790 MeV/nucleon, plotted against the occupation 


probability of the valence neutrons configuration in p, d and s orbit. 


TABLE 3. Two-neutron occupancy (96) for 12Be (g.s.) as 10Be + 2n, 


(1p)? (1d)? (25)? Method Ref. 
38 29 32 Calc. [20] 
74 18 8 Calc. [51] 
32 15 53 Calc. [21] 
13-19 10-13 67-76 Calc. [24] 
31 34 35 Calc. [22] 
25 (3He,5He) [52] 
20 B decay [53] 
35 B decay [54] 
50 8 decay [55] 
14-20 nucleon transfer reaction [3] 
32 30 38 extracted from spectroscopic factors [4] 
25 «5 extracted from Gamow-Teller Transition Strengths [10] 
« 5696 extracted from o present work 


V. SUMMARY 


In summary, the ground state structure of !?Be is of great 
significance. The ambiguity of the configuration of !?Be 
ground state has attracted our attention. We bring forward 
to determine the s orbital component of !?Be ground state 
by measuring the c; of !?Be on C and AI at several tens of 
MeV/nucleon. By using existed interaction cross section data 
of !?Be on C at 790 MeV/nucleon, we roughly determine the 
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